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Abstract 

This paper describes a novel directional neutron detector prototype. The low 
pressure time projection chamber uses a mix of helium and CF4 gases. The 
detector reconstructs the energy and angular distribution of fast neutron recoils. 
This paper reports results of energy calibration using an a source and angular 
reconstruction studies using a collimated neutron source. The best performance 
is obtained with a 12.5% CF4 gas mixture. At low energies the target for fast 
neutrons transitions is primarily helium, while at higher energies, the fluorine 
contributes as a target. The reconstruction efficiency is both energy and target 
dependent. For neutrons with energies less than 20 MeV, the reconstruction 
efficiency is ~40% for fluorine recoils and ~60% for helium recoils. 



1. Introduction 

Precise measurement of the energy and direction of cosmogenic-muon-induced 
fast neutrons, with various levels of underground shielding, will be valuable for 
future neutrino and dark matter experiments pQ. In this paper, we present a 
prototype detector for this purpose. The detector was inspired by the design for 
the directional dark matter search DMTPC [2]. We introduce modest modifica- 
tions for directional neutron detection [3]. Our first application of this detector 
will be in the Double Chooz (DC) neutrino experiment [4], and it is, therefore, 
called DCTPC. 

DCTPC is a low pressure time projection chamber (TPC) filled with two 
gasses: (1) He, which is the fast neutron target, and (2) CF4 which serves as 
a scintillant and a quencher for the electron avalanche. He and F have similar 
neutron cross sections, but the He nucleus acquires more energy and leaves a 
longer trail of ionization along its track. This makes the He recoil easier to 
reconstruct. The directionality of the recoil can be used to discriminate against 
non-cosmogenic neutron backgrounds and isolate the signal. This detector is 
blind to X-rays and minimum ionizing particles (MIP), like muons, because the 
density of primary ionization electrons from MIP particles is too small to be 
detected. 

In the DCTPC design, the ionized electrons from a particle track drift down 
to a stainless steel ground mesh. Between the ground mesh and the anode, an 
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Figure 1: The recoil energy spectra of neutrons from the KARMEN2 shielding \E\ (top) and 
the rock surrounding LVD 6 (bottom). These distributions include the detector response. 
Correcting for the kinematic difference between neutron recoils in He gas and recoils in liquid 
scintillator, DCTPC can probe the three populations neutrons that create recoils of approxi- 
mately <20 MeV in these experiments. 



avalanche of electrons occurs. A charge-coupled device (CCD) camera installed 
in the detector images the visible scintillation light from the avalanche. The 
charge created in the avalanche is readout from the anode plate and ground 
mesh. The track direction and energy can be reconstructed using the CCD 
image and the charge information. The reconstructed track from the recoiling 
nucleus can then be correlated to the incoming neutron direction and energy. 

In this article we present on-surface measurements from a small 2.8 L proto- 
type of DCTPC. We provide results from an a-source run and a 252 Cf neutron 
source run as well as an on-surface cosmogenic run. These results demonstrate 
the capability of the design for performing an underground measurement. These 
source- studies, coupled with a future underground measurement, will guide the 
design for a larger multi-detector system for Double Chooz. 



2. Motivation for DCTPC 

DCTPC is under construction for the Double Chooz neutrino oscillation ex- 
periment [4] located at the Chooz Nuclear Power Station, a site with two 4.27 
GW thermal reactors. Double Chooz has a near-far design with distances to 
reactor 1 (2) of 456.5 m (339.2 m) for the near site and 1114.7 (998.0 m) for the 
far site. The detectors sites have 114 m.w.e and 300 m.w.e shielding, respec- 
tively. Antineutrinos from the reactors are detected in Gd-doped scintillator oil 
targets. Fast neutron interactions are a background to the antineutrino events 
used in the oscillation search, and so must be well-understood. 
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DCTPC provides a benchmark for the Double Chooz fast neutron simulation 
at the specific sites of the two Double Chooz detectors. The neutron flux energy 
and angle measurements from the rock can be used to tune the Monte Carlo 
and demonstrate a clear understanding of this background in the Double Chooz 
neutrino oscillation analysis. 

However, the motivation for DCTPC goes beyond its specific applications in 
Double Chooz. Backgrounds from fast neutrons produced in rock are an issue 
for all low background experiments. Below, we review the expected sources of 
fast neutrons and the existing measurements. We point out where DCTPC can 
provide unique data. This includes energy spectra and angular distribution as 
well as measurement of total rates. 

The final DCTPC design calls for a large detector located at each site. The 
small prototype described here will be installed in the 300 m.w.e site for ini- 
tial studies. These shallow depths are particularly interesting for studies of 
environmental effects on neutron rates that are relevant for other shallow-site 
experiments. 

2.1. Measuring Three Sources of Neutrons at Depth 

Fast neutron sources from the rock environment of underground experiments 
fall into three categories summarized in Fig. [I] On this figure, we also show the 
results from experiments at two depths, KARMEN2 [5 (surface) and LVD [6] 
(3200 m.w.e.). 

The lowest energy category arises from neutrons produced by (a, n) interac- 
tions. The a's are produced by the natural radioactivity in the rock surround- 
ing the experiment. At locations with low shielding, this neutron source is 
overwhelmed by neutron production from cosmic ray muons. However, muons 
are attenuated with depth and the (a, n) background becomes significant at 
300 m.w.e. [7]. DCTPC, with its very low energy threshold (80% charge re- 
construction efficiency at 200 keVee, as shown below), can sample these (a, n) 
events. The comparison of the 114 m.w.e. to 300 m.w.e. depths at Double 
Chooz will be useful in isolating the (a, n) source from the muon-induced neu- 
trons. The (a, n) measurements will be unique to the U/Th content of the 
rock surrounding Double Chooz, however they will be useful for benchmarking 
simulations for wider use. 

The neutron sources from cosmic ray muon interactions from rock fall into 
two broad categories: stopped muon capture and deep inelastic interactions. 
The rates depend upon the depth of the experiment because low energy muons 
are attenuated with shielding. DCTPC sensitivity will reach to >20 MeV in fast 
neutron reconstruction, overlapping with the KARMEN2 and LVD data sets. 
As we show below, the maximum reconstructed recoil energy depends upon the 
size of the detector. 

2.2. Total-Rate Measurements for Benchmarking Simulations 

Most experiments, including Double Chooz, rely on simulations, typically 
FLUKA [3 [9] and GEANT4 [TOl fTT]. to determine the expected fast neutron 
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rate. These simulations generally focus on the > 10 MeV fast neutrons that are 
produced by muons. The uncertainties on these calculations are large and more 
data will be useful as a constraint. 

The exercises in benchmarking simulations [TJ [12] [13j [7J [14j [15] mainly em- 
ploy data on neutrons produced by muons passing through scintillator-detectors 
installed at various depths. This is a relatively well-controlled circumstance for 
data-to-monte carlo comparison, since one can trigger on the muon and then 
study the subsequent neutron interactions. However, in this case the muon- 
targets are carbon and hydrogen, rather than rock. This results in uncertainty 
from nuclear effects when the interaction rates are translated to the background 
rates from the heavier isotopes that make up the rock at a particular experi- 
mental site. 

Measurements of the total rates at various depths on hydrocarbons include 
Refs. [6j [16] [TTJ [THJ [19] [20] [21] . In principle, neutron production by muons in hy- 
drocarbons should be easier to simulate than neutrons produced in a surround- 
ing, conglomerate-material rock-environment, as discussed above. However, in 
practice, agreement between data and simulation for oil-produced neutrons is 
poor [12j [13] [TJ [TJ [15] . DCTPC adds data points at two depths for benchmark- 
ing the total rate directly from the rock, which will complement the existing 
measurements on scintillator. 

2.3. Energy and Angle Dependent Data 

Compared to the total-rate measurements, information on the energy spectra 
of fast neutrons is scarce. Data on the energy spectra have been published by 
only three experiments: These are the KARMEN2 experiment on the surface [5 , 
KamLAND at 2700 m.w.e. [20^ and the LVD experiment at 3200 m.w.e. [6]. 

The KARMEN2 and LVD spectra are shown in Fig.[l] Both KARMEN2 and 
LVD used a muon trigger, and so neither data set extends into the low energy 
(a, n) interaction region. DCTPC will nicely fill the gap in data as a function of 
depth between KARMEN2 and the deeper experiments. There are important 
differences between the KARMEN2 and LVD data sets, aside from depth, that 
must be considered when cross-comparing and using these as benchmarks. The 
LVD neutron spectrum is dominated by neutrons from muon interactions on 
the scintillator. Therefore extrapolations back to the spectrum from rock face 
have the same uncertainties as scaling the total production rates. The KAR- 
MEN2 spectrum is from neutrons produced in the steel shielding surrounding 
their detector. The KARMEN2 shielding was instrumented with thin scintilla- 
tor counters, which provided the muon trigger for this study. This steel target 
is more similar to rock environments underground, making this data attrac- 
tive for study. However, the fact that this is an iron rather than hydrocarbon 
target is not accounted for in the benchmark from Ref. [15 . The KamLAND 
measurement is also not straightforward to include in comparisons. 

The KamLAND data set is significantly different from KARMEN2 and LVD, 
and so is not shown in Fig. [I] KamLAND isolated events by requiring a fast- 
neutron-like scatter in an outer detector and a neutron-capture-signal in an 
inner detector region. This analysis measures incoming neutrons from rock, 
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which is ideal. However, the neutrons must penetrate substantial outer regions 
of the detector, including a water veto and stainless steel sphere, in order to 
reach the inner fiducial volume for measurement. In the process, the neutrons 
are likely to scatter more than once. As a result, the energy spectrum of the 
neutrons is substantially modified, and interpretation of the data requires a full 
detector Monte Carlo simulation of KamLAND. 

One might argue that simulations are best tuned from accelerator-based 
experiments, where one has more control over the target material and muon 
energies. Ref. [14] analyzed the results of CERN NA55 [22 and a SLAC beam 
dump [23j [24] experiment. They conclude that, even within this more con- 
trolled environment discrepancies exist between the data and the simulation. 
Re- weighting the simulation corrects normalization discrepancies. However dis- 
agreements with the energy distribution remain. 

Along with a measurement of the energy spectrum, DCTPC also provides 
an angular distribution of events. In this paper, we show the quality of this 
angular reconstruction. The unique directional information can help disentangle 
different neutron sources. 

2.4- Studying Environmental Effects at Low Shielding Depths 

Double Chooz is located at relatively low shielding depths compared to most 
underground experiments. However, there are a number of projects running or 
planned at modest depths. These include COUPP [25 and other non-accelerator 
experiments that run in the NuMI near hall (^300 mwe) at Fermilab. The 
proposed Homestake LAr detector currently assumes a hall below only 800 ft 
of rock [26]. There is also a Japanese proposal to build a KamLAND-style 
geoneutrino detector most likely located at 300 m.w.e. [27 . Therefore, there 
is motivation to understand environmental effects that may be unique to labs 
with modest overburden. 

In particular, there is some evidence that the fast neutron flux varies with 
rainfall. The increased water level acts as a moderator for both the cosmogenic 
muons and their daughter neutrons. Fig. [2] shows fast neutron and rainfall data 
from the Baradello Underground Laboratory. The lab is located in Italy at a 
depth of 300 m.w.e. [28]. If the rainfall data is offset by 30 days, reductions in 
the neutron rate track increases in rainfall well and vice versa. In addition to the 
rate, presumably the energy and direction of the fast neutrons are also affected 
by the water, although the Baradello detector cannot provide this information. 
It would be valuable to repeat this study, and eventually to understand if this 
effect persists at greater depths. 

Other laboratories, including Double Chooz, are likely to see a correlation 
between rainfall and neutron rates, although with a different relative time offset 
because of differences in rock structure and water flow. This motivates develop- 
ment of a detector which measures the energy, direction and rate of fast neutrons 
produced in the Double Chooz detector hall, as a function of time. This can be 
used to study time- variations in the neutron flux with the physical environment 
and accurately tune the rate, energy, and direction of fast neutron simulations 
including time dependent effects. 
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Figure 2: Neutron rates per day (black symbols), and rainfall per day shifted by 30 days (red 
line) at the Baradello Lab located at 300 m.w.e. 



3. The Prototype Detector 

We have developed a prototype of DCTPC, and in this section we describe 
the design of this detector. The prototype is a TPC with a drift volume above 
an amplification plane that has two active regions. The amplification plane 
is imaged with one CCD, and the charge is read out for each active region. 
The redundancy in the readout of the amplification plane is powerful for event 
reconstruction and background rejection. The detector schematic is shown in 
the Fig. [3] and a photograph of the detector vessel is shown in Fig. |1J 

3.1. The Time Projection Chamber 

The detector follows the classic design of TPCs. We have chosen a mixture of 
CF4 and He as our operating gas. We study several gas ratios below, but most 
data were taken with the gasses mixed at 75 Torr and 525 Torr respectively. 
Helium and fluorine both have large cross sections for fast neutron scattering 
[29] . The 4 He nucleus is closer to the neutron mass, and therefore the recoils will 
produce longer tracks than the fluorine nucleus. Arguably, H2 is a better match 
to the neutron mass and, in that sense, is a better target. However, we prefer 
to use inert gasses which allow operation without substantial safety concerns. 
In addition H2 has a lower cross section and produces less ionization along the 
track. 

Within the vessel, a uniform drift field is molded by a series of copper rings 
with an inner diameter of 26.9 cm. The rings are connected with resistors, see 
Fig. [3j The top ring carries a high transparency stainless-steel mesh with wire 
pitch of 512 /am and wire diameter of 31 jam. This top mesh is maintained at 
a negative high voltage. The total height of the field cage is 10 cm. Primary 
ionization electrons are drifted down to a grounded stainless-steel mesh with 
wire pitch of 256 /am and wire diameter of 30 jam. 
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Figure 3: A schematic of the detector: the drift field is created by a cathode mesh, field- 
shaping rings attached to a resistor chain, and a ground mesh. Primary ionization from a 
recoiling nucleus is drifted down to the ground mesh. The high-field amplification region is 
formed by the ground mesh and the anode plane. The ground mesh is read out with a fast 
amplifier and the veto and anode are readout with charge-sensitive preamplifiers. Scintillation 
light from the amplification region is recorded with the CCD camera. 
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Figure 4: Left: The prototype detector. Left: The detector vessel is sitting on top of the 
cart which holds the readout electronics. The CCD camera is mounted on the top flange. 
Right-top: A diagram of the anode plane. The central anode region is surrounded by the 
veto ring, which is then surrounded by a grounded ring. Spacers sit on top of the anode 
plane to preserve the space between the anode plane and the ground mesh. Right-bottom: A 
photograph of the field cage and amplification plane of the TPC. The ground mesh is difficult 
to see, but the spacers separating it from the anode plane are visible as horizontal lines. 
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The amplification region is created by the ground mesh and the copper 
anode that are separated with insulating tubes 440 fim in diameter. These 
insulating tubes are called spacers and are placed 2.56 cm apart [30]. This high- 
field amplification gap multiplies the primary ionization charge by a factor of 
approximately 10 5 , as we discuss below. Operating voltages are gas-dependent. 
The anode of the amplification plane is divided into two regions. The inner 
part of the anode, diameter 24.7 cm, is used for the detection of tracks. A veto 
ring, with inner diameter 24.8 cm and outer diameter 26.8 cm, surrounds the 
central plane and is used to veto charged particles originating from decays from 
within the detector components. A central square region, 16.7 cm x 16.7 cm, of 
the amplification plane is read out optically using the CCD camera. Using this 
geometry, we calculate that the active volume of the detector that has a charge 
readout is 5.6 L, while the volume with both CCD and charge readout is 2.8 L. 

3.2. CCD Readout 

The CCD camera records scintillation light created in the amplification re- 
gion, see Fig. [3] The camera is the U6 model made by Apogee Instruments Inc., 
with the Kodak KAF-1001E CCD chip. The chip is kept at -20 °C to minimize 
the dark current. An image from the amplification plane is focused to the CCD 
chip using a Nikon lens, of focal length 55mm and set to f/1.2. The CCD chip 
consists of 1024 by 1024 square pixels where the side length of a pixel is 24 jira. 
The active area on the anode plate is a square of 16.7 cm x 16.7 cm, see Fig.[4j 
In order to maximize the signal-to-noise ratio, we merge the CCD pixels into 
4x4 bins during readout, so that the effective spatial resolution is 640 fim. 

The camera does not use a shutter and the CCD sensor is always live, even 
during the CCD readout. The unit for one pixel count returned by the CCD 
readout is an ADU (Analog-to-Digital Unit). The fact that the sensor is always 
live creates two types of exposures: a nominal exposure lasting 1 second, and 
a 'parasitic' exposure occurring during the CCD readout and event processing. 
The 'parasitic' exposure lasts approximately 250 ms. Track images recorded 
during this time can be shifted in the direction of readout of CCD rows, and do 
not have corresponding signal in the charge readout channels. Therefore, they 
are easy to remove during analysis. 

3.3. Charge Readout 

The anode plate contains three regions separated by electrically insulated 
gaps, as shown in Fig|4] The outermost circular region is connected to a common 
ground, onto which the ground mesh is attached. The inner disk region is the 
actual anode, maintained at 680 V. The circular-ring region between is used as 
a veto and is kept at the same potential as the anode. 

When a particle enters from outside, its path must cross the veto region. 
Therefore, we can detect this type of event by inspecting the charge signal from 
the veto region. The veto is not applied at the trigger level. The veto signal is 
stored for use in the analysis. 

The anode signal is collected using a charge-sensitive preamplifier (CSP), 
the Cremat CR-113 with a gain of 1.5 mV/pC. The pulse height of the output 
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of the CSP corresponds to the total charge collected on the anode. Therefore, 
the pulse height is used to reconstruct the total energy of the event. The veto 
electrode is read out through the more sensitive CSP, the Cremat CR-112 with a 
gain of 15 mV/pC. The ground mesh is read out through a fast current amplifier, 
the Route2Electronics HS-AMP-CF-2nF with 1 ns rise time and gain of 80. The 
fast amplifier has a built-in protection that guards against discharges from this 
high-capacitance detector. The CSPs are protected with 300 ft resistors. All 
three signals are digitized using the ATS860, a 250 MS/s 8-bit flash ADC made 
by AlazarTech. The digitization of all three channels is triggered by a signal on 
the anode and stored with the current CCD image. The trigger threshold used 
is gas dependent. 

The charge readout is very powerful in suppressing CCD-specific backgrounds 
(e.g., direct hits in CCD chip and residual bulk images). Introduction of the 
charge readout is the major improvement over the CCD-only approach described 
in Ref. [3]. 

3.4- Event Readout, Reconstruction and Correlation 

As discussed above, the CCD camera is constantly taking data in sequential 
exposures of 1 s with effectively zero dead-time. We call this mode of data 
taking "sequential read-out" in the discussion below. The rate of events in 
the detector, even with the neutron and a sources, is relatively low, and most 
exposures contain no tracks. Since tracks imaged by the camera do not have 
an exact time stamp, those CCD images with tracks are matched to the proper 
charge-readout sub-event using the reconstructed energy. In this way, the track 
is assigned a timestamp. Because charge-readout has dead-time due to the 
digitization, there are some fraction of images that are missing the corresponding 
charge event. 

The combination of the four measurements, CCD image, charge deposited 
on the inner anode and veto anode and the current recorded at the ground mesh, 
is used for track reconstruction. Given the value of the electron drift velocity 
in a particular gas, the pulse rise time of a waveform from the ground mesh can 
be used to reconstruct the length of the track in the drift direction, z for this 
detector. The CCD image provides a 2-dimensional projection of the track on 
a horizontal surface, x-y. We can then use these two pieces of information to 
reconstruct the 3-dimensional length of a track in each event. Unfortunately, the 
drift voltage is not sufficient, and diffusion dominates the rise time of pulses from 
the ground mesh. For the following analyses, we will only use the 2-dimensional 
tracks obtained from the CCD image. 

The CCD track reconstruction starts with subtracting the bias image that is 
obtained at the beginning of every run. At this time, any unusually high-valued 
pixels are removed. The image is then smoothed. With the smoothed image, a 
search is conducted for clusters of adjacent pixels with values 3.7a larger than 
the mean value of the image [2 . If enough such pixels are found, the cluster 
is determined to be a possible track. An additional ring of pixels around the 
main cluster is added to account for diffusion, and the track is the corresponding 
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pixels from the original un-smoothed image. For each of these tracks, the energy, 
range, angle, maximum value, and several moments are calculated. 

In the following analyses, we require that the track be contained within the 
TPC. The veto allows us to remove background radiation from the detector 
itself, for example a particles emitted from the drift rings or the vessel wall. 
However, this also vetoes long tracks produced within the active volume which 
happen to enter the veto. Additional inefficiency in the reconstruction of long 
tracks comes from the requirement that the track be fully contained in the CCD 
view- field. We show below that this is not a significant constraint, and we discuss 
the reach in neutron energies that we can study with the small prototype. This 
will be even less of an issue in the final DCTPC design, which will be larger, 
and allow for longer tracks. 

4. Spatial Calibration 

The spatial calibration is obtained by measuring the positions of the wire 
spacers separating the anode plane from the ground mesh, shown in Fig. |4j The 
spacers create a region of low gain. A 137 Cs source is used to generate electronic 
recoils so that a small but diffuse amount of light is seen in each image. By 
integrating many images, the spacers become visible as regions with low light 
levels. The positions are determined by projecting the integrated image onto 
the axis perpendicular to the spacers, so that the spacers appear as Gaussian 
peaks superimposed on a background. The mean positions of the peaks are fit 
to a line, and the slope gives the separation between peaks. The spacers are 
placed 2.5 cm apart, and the fit gives a calibration of 163 fim per pixel. The 
total area of the 1024 by 1024 pixel region imaged by the CCD is then 16.7 cm 
x 16.7 cm, as discussed in Sec. |3.2| The uncertainty on the spatial resolution 
is believed to be dominated by systematics from any misplacement of spacers, 
optical distortions, and deviations in the spacer shape from a line. 

5. Energy Calibration 

The energy calibration is performed with a 241 Am a source. The source is 
contained in a tube which collimates the a particles. This insures that the a 
tracks are measured by both the CCD and the charge readout. A window on 
the source degrades the energy of the a particles leading to a mean energy of 
4.40 ± 0.04 MeV [3T] , Example a events in three different gas mixtures are 
shown in Fig. [5] 

The total ionization left by the track is measured in three ways: (1) inte- 
grating the light in the CCD track, (2) measuring the pulse peak height of the 
central anode signal and (3) integrating the current signal from the mesh chan- 
nel. At the drift voltage used in these data, the a tracks saturate the digitizer 
of the ground mesh preamp and possibly the preamp as well. For this reason 
only methods (1) and (2) are used for determining track energy. The raw units 
of CCD energy are ADU and those for the charge readout are mV. 
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Figure 5: Examples of a particle tracks from the 241 Am source. From the top to bottom, 
graphical images of the CCD and charge read-out system are shown. Top: CCD image of 
a particle track. Middle: Waveforms from the Veto and Anode charge read-out. Bottom: 
Waveforms from the Ground-Mesh charge read-out. Each column corresponds to a different 
gas mixture. Left: 100% CF4. Middle: 6.25% CF4 + Helium. Right: 12.5% CF4 + Helium 



The energy distributions from the anode and CCD are shown for all three 
gas mixtures in Fig. [6j The distributions of track CCD energies are fit to a 
Gaussian distribution, with the mean value giving the energy calibration. The 
energy distributions for the charge signals shows some distortion, and these are 
best fit using a log-normal distribution. The results of the calibration are given 
in Table [T] Since quenching of the a track is a negligible ~1% effect, the energy 
is quoted in units of electron equivalent energy keVee without any additional 
corrections. 

We can validate our calibration method by comparing the correlation be- 
tween the CCD energy and the summed energy of anode waveforms during the 
corresponding CCD image. As discussed in Sec. |3.4[ the sequential CCD camera 
read-out produces images that may contain multiple tracks or, in most cases, no 
tracks at all. We can confirm a linear relationship between the summed anode 
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Figure 6: Energy gain calibration using the charge read-out (top) and CCD (bottom). From 
the left, 100%, 12.5%, and 6.25% CF4 mixtures with He are shown. 
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Table 1: Summary of detector calibrations for running in the three different gases. The trigger 
threshold is that used for neutron and background running. Higher thresholds are needed for 
the energy calibration runs. 



energy and the CCD energy. This implies that a good match between a track 
and an anode waveform can be found for most of events. This is shown in Fig. [7] 
The gas gain is calculated from the charge energy calibration, the preampli- 
fier gain of 1.5 mV / pC, and the gas work function. The work function of CF4 
gas was measured to be 33.8 eV, see Ref. [32]. This is sufficient to determine 
the gain in pure CF4. The work functions for the helium gas mixtures are not 
known. In the gas mixtures, the nuclear recoils will primarily ionize the CF4, 
so this value is a good approximation. Helium is known to have a larger work 
function than CF4, so the true gain in the gas mixtures, summarized in Table [I] 
will be slightly larger than the value reported here. 



6. Selection Cuts 

In order to identify the tracks from nuclear recoils with energies much lower 
than those from the a source's tracks, a number of selection cuts are applied to 
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Figure 7: Comparison of energy from CCD light versus charge read-out signal. Left: 100% 
CF4. Right: 12.5% CF4 + He. 



remove various backgrounds. These backgrounds include CCD noise artifacts, 
electronic recoil charge signals, and ionization from particles directly passing 
through the CCD. Different cuts must be applied to each gas in order to properly 
filter out these background events. 

6.1. CCD Cuts 

The CCD cuts for each gas are very similar. All CCD tracks are required 
to have greater than 1 ADU. The tracks are required to have a nonzero range, 
meaning at least two adjacent pixels have counts at least 3.7 a above the mean 
for the image [2 . Tracks must not reach the edge of the image, and in pixel 
coordinates this means 24 < x, y < 1000. This ensures that the track is fully 
contained within the image. Ionization within the CCD chip is removed by 
requiring the maximum pixel value in a track to be less than 250 ADU. Noise 
artifacts are removed by requiring that the maximum pixel in the image contains 
less than 25% of the total signal. Hot pixels and residual bulk images are rejected 
by requiring that no two tracks occur within 12 pixels within a run. The event 
rate, when no sources are used, is low enough that few real tracks are rejected by 
this cut. Finally, a gas dependent maximum range is used to reject background 
a events. In CF4, tracks must be shorter than 80 pixels (1.3 cm). In the 12.5% 
mixture and 6.25% mixture, the values are 160 pixels and 140 pixels respectively. 
The maximum range in the 6.25% mixture is limited by the higher gain and the 
maximum energy that can be measured before saturation becomes an issue. 

6.2. Charge Cuts 

The charge cuts are more gas dependent than the CCD cuts. The drift 
voltage in pure CF4 is closer to the ideal value for minimizing diffusion. This 
creates good separation between the peaks from the electron avalanche and the 
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% of CF 4 


tyeto 


tanode and ^mes/i 


Pulseheight Range 


100% 


<800 ns 


tmesh < 30 ns for pulseheight < 100 mV 
tmesh < 25 ns for pulseheight > 100 mV 


25% to 75% 
25% to 75% 


12.5% 


<800 ns 


tanode < 1250 ns 


10% to 90% 


6.25% 


<900 ns 


tanode < 740 ns 
tanode < 1400 ns 

tmesh < 600 ns 
tmesh < 800 ns 


25% to 75% 
10% to 90% 
25% to 75% 
10% to 90% 



Table 2: Rise time requirements. Variables are defined in the text. Events are accepted if they 
pass the requirements listed. The pulseheight range over which t ve to is measured is always 
10% to 90%. The pulseheight range for the definition of tmesh an d t ano d e vary as listed in 
column 4. 



% of CF 4 


Ryeto 


Relec 


Rion 


Rmesh 


100% 


>4 


1.7 < Relec < 2.7 


1.9 < R ion < 3.5 


N/A 


12.5% 


>3 


N/A 


N/A 


1.3 <C Rmesh ^ 3 


6.25% 


>4 


N/A 


N/A 


N/A 



Table 3: Pulse height ratio requirements. The ratios are defined in the text. Events are 
accepted if they pass the requirements listed in this table. N/A means that a reliable cut on 
this variable could not be identified. 



drifting ions in the ground mesh pulses. In the helium gas mixtures, the drift 
field is much less than the ideal value, so much more diffusion is seen. Separate 
electron and ion peaks within a pulse are rarely seen, thus only the overall pulse 
height and rise times are used in the analysis of the veto, anode, and ground 
mesh signals. The significant differences in the responses of the three gasses is 
reflected in the cuts used to retain events, Tables [2] and [3j 

The charge analysis uses several variables related to the pulse shape. For 
the veto we define two variables: the time is takes the pulse to rise from 10% 
to 90% of the peak pulse height, t veto , and the ratio of the mesh peak pulse 
height to the veto peak pulse height, R ve t . For the mesh and anode, we define 
pulse rise times, tmesh an d ^anode- The pulse height range over which tmesh an d 
tanode are defined is tuned for each gas mixture. For the mesh we also define 
two ratios for the pulse height: R e iec the ratio of the electron peak mesh pulse 
height to the height of the anode pulse and Ri on the ratio of the ion peak mesh 
pulse height to the anode pulse height. Rmesh is a simple ratio of the anode 
pulse height to the mesh pulse height. 

These cuts are tuned to reject background due to electronic recoils and cross 
talk between the anode and the veto. The tuning of the cuts to retain signal and 
reject these backgrounds highlights the difficulties with using helium and the 
accompanying higher background rates. These difficulties are balanced by the 
advantages of longer recoil tracks from this target nucleus. The gas mixtures 
with helium have a substantially slower rise time for the signal and this is 
reflected in the requirements of Table [2] For the pure CF4 gas mixture, the 
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electron and ion signals can be separated, allowing cuts on R e \ ec and Rion, as 
listed in Table [3| However, in the 12.5% mixture we cannot identify reliable 
cuts for these variables and must use R m esh to isolate the signal. In the 6.25% 
mixture, we could not identify any pulse height ratio cuts that efficiently selected 
signal, and therefore additional rise time cuts are used. 

6.3. Charge-Light Matching 

In the data presented here, only events passing all cuts described in Sec- 
tion |6.1| and containing a single CCD track are used. In normal running it is 
unlikely to see events with multiple true nuclear recoils, so this later cut has a 
negligible effect on the efficiency. The energies of the charge signals which pass 
the cuts described in Section [6^2] are compared to the energy of the CCD track, 
and the signal with the best energy match is identified as the true signal from 
the track. The energy from the anode channel is expected to be more accurate 
than the CCD energy, so it is used as the final track energy. 

In pure CF4, the anode and CCD energies must differ by less than 40 keVee 
if the anode energy is less than 1 keVee and by less than 70 keVee otherwise. 
In 12.5% CF 4 , the difference must be less than 75 keVee for E ano ^ e < 250 
keVee, 125 keVee for 250 < E ano d e < 500 keVee, and less than 150 keVee for 
E anode > 500 keVee. We did not set charge matching cuts for the 6.25% mix. In 
this gas, the absolute and relative energy calibrations of the different channels 
are found to be inconsistent over several gas fills, likely due to the much steeper 
increase of gain with respect to anode voltage. In this gas, small variations in 
voltage can cause very large gain differences, and it was determined that this 
gas is not to be stable enough to be useful in the current detector configuration. 
The efficiency of the charge cuts for the other two gases, including charge-light 
matching is given in Fig. [9j This is done with neutron source data described 
below. This data is dominated by nuclear recoils rather than CCD artifacts, so is 
ideal for the efficiency determination. The efficiency is defined as the ratio of the 
number of tracks passing all CCD cuts, charge cuts and charge-light matching 
compared to the number of CCD tracks found without using the charge signals. 

7. Neutron Source Run 

The response of the detector to neutrons was characterized using a 1.24 mCi 
252 Cf source. 252 Cf decays via spontaneous fission 3% of the time emitting 
multiple neutrons with a mean energy of 2.35 Me V [33 . Given the activity of 
this source, we expect 5.7 x 10 6 neutrons per second in all directions. The source 
is placed in a collimator made of borated plastic, 2.1 m from the detector. The 
high flux of both neutrons and x-rays causes an increase in the rate of sparks 
in the TPC from ~0.01 Hz to 0.1 Hz. A spark is identified as an image with 
a mean number of ADU 1% greater than the previous image. Events within 6 
exposures of a spark are excluded from analysis to allow sufficient time for the 
high voltage to recover. 

The neutron source generates a large number of nuclear recoils, with 0(1) 
neutron-induced nuclear recoils detected per 1 s exposure. A large number of 
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Figure 8: Example figures of neutron recoils from 252 Cf source. From the top to bottom, 
graphical images of CCD and charge read-out system are shown. Top: CCD image of 4 He 
recoil. Middle: Waveform from the Veto and Anode charge read-out. Bottom: Waveform 
from the Ground-Mesh charge read-out. Each column corresponds to a different gas mixture. 
Left: 100% CF4. Middle: 6.25% CF4 + Helium. Right: 12.5% CF4 + Helium. 
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Figure 9: Efficiency of charge analysis in accepting tracks identified by the CCD analysis. 
The efficiency plateaus at 70-80% are due to exposure difference between the charge and light 
readout and depend on the event readout time. The decrease in efficiency at higher energies is 
from the mesh pulse reaching the maximum value read out by the digitizer. Left: 100% CF4 
gas. Right: 12.5% CF4 + He mixture. The errors on the data points reflect the statistics of 
the calibration run. 



electronic recoils from the high 7 flux are observed in the charge signals. How- 
ever the CCD remains blind to these events due to their low primary ionization 
density. Additionally, the source creates a much higher number of CCD-induced 
artifacts, likely from electronic recoils inside the active volume of the CCD chip. 
Even with these large source activity issues, the cuts described in section Sec- 
tion [6]T] and Section 6.2 are sufficient for this analysis. 



After applying all cuts, we use the measured track angles and ranges, Fig. 10 



and Fig. [TT] to evaluate the performance of the detector and the analysis. The 
measured two-dimensional track range generally falls just below the expected 
mean three-dimensional range calculated with the publicly available program 
SRIM [34 . For the 12.5% CF4 mixture, the measured two-dimensional track 
range falls close to the expected range for helium tracks. This follows from 
the fact that helium recoils are kinematically favored across most of the energy 
range measured here. At lower energies, where both helium and fluorine recoils 
are expected, the track ranges appear to be consistent with either nucleus. 

We use the track angles to reconstruct the location of the neutron source. 
Tracks along the direction of the incoming neutrons should have a decreasing 
light profile, and this "head-tail effect" can be used to reconstruct the track 



direction [35]. Two peaks are seen in the angular distribution of Fig. 10 The 
larger peak at —90 degrees is along the true mean direction of the nuclear 
recoils. The second peak at +90 degrees is primarily from events where the 
track direction was mis-reconstructed. The large width of the distribution is 
expected from elastic scattering kinematics. Additionally, some nuclear recoils 
are from neutrons that have scattered several times within the laboratory, and 
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no longer appear to originate at the source position. 

A model of helium and fluorine recoils is fitted to the measured recoil spec- 



trum in Fig. 12 The fit uses the likelihood approach with the probability 
density functions for helium and fluorine based on simulation, and a 252 Cf spec- 
trum based on ENDF tables [29]. We find the fraction of helium recoil events 
to be 61 ± 1 %. Another way to present this result is in terms of the unfolded 
neutron spectrum. We use the recoil data and the unfolding matrix, based on 



the simulation, to obtain the neutron spectrum as shown in Fig. 12 There is 



a reasonable agreement with the neutron spectrum from a 252 Cf source that is 
based on ENDF tables [29 . In order to improve the agreement between the data 
and the simulation of the recoil and neutron spectra, we need to improve the 
simulation to include the efficiency loss due to sparks. This requires additional 
measurements, and a full simulation of the detector material which modifies the 
energy spectrum of neutrons entering the sensitive region. 



8. Surface Run 

With the DCTPC prototype, we have also taken data with no sources to 
measure cosmogenic, environmental and detector backgrounds in a surface lab 
on the MIT campus. The lab is on the ground floor of a two-story building. 
The lab has approximately 1 m concrete walls on all sides. After removing 
spark events, we obtain 112713 exposures in pure CF4 and 103822 exposures 
in the 12.5% CF4 mixture. As in the other runs, each exposure is 1 s, so this 
corresponds to 1.3 and 1.2 days of exposure, respectively. 

In the pure CF4 run, we identified 128 possible tracks in the CCD analysis. 
Most of these tracks appear to be from CCD backgrounds such as noise artifacts 
or ionization in the CCD. Of these 128 tracks, only 10 were matched to a 
charge trigger. Several additional tracks had too much energy to be accurately 
measured by the charge channels, although the pulses were collected. These 
tracks are mostly low energy [E < 100 keVee) and the ranges generally appear 
to be consistent with the expected values from SRIM for fluorine recoils. A 
single event has a range more consistent with helium or hydrogen tracks. 

In the 12.5% mixture run, we identified 71 potential tracks in the CCD anal- 
ysis, with 16 having a matching charge signal. As with the pure CF4 run, most 
of the potential tracks seen in the CCD are likely noise artifacts or tracks left 
directly in the silicon of the CCD chip. As expected, we see more tracks in this 
run due to the presence of helium. The tracks seen here are also generally at 
higher energies than were seen in pure CF4. This is to be expected if most of 
the tracks come from elastic scattering from neutrons due to the favorable kine- 
matics of neutron-helium scattering compared to neutron-fluorine scattering. 
Fewer events at lower energies are seen due to the differing energy thresholds of 
the two gas mixtures. There does not appear to be a favored direction of these 
tracks, although much higher statistics would be necessary to confirm this. The 
events with E > 300 keVee here have ranges consistent with helium recoils but 
not carbon or fluorine, while most of the other tracks have ranges consistent 
with either carbon or fluorine nuclei or near vertical helium nuclei. 
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Figure 10: Energy-Angular distribution for calibration and background runs different gas 
mixtures. Top: Calibration run with 252 Cf source. Bottom: Background run. Left Column: 
100% CF4 gas. Right Column: 12.5% CF4 + He mixture. The direction of the neutrons from 
the source is indicated by a red solid line at —90 degrees. The population at +90 degrees is 
due to incorrect reconstruction of the track direction. 
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Figure 11: Energy- Range distribution, calculated in two dimensions, for calibration and 
background runs with two different gas mixtures. Top: Calibration run with 252 Cf source. 
Bottom: Background run. Left Column: 100% CF4 gas. Right Column: 12.5% CF4 + He 
mixture. Three-dimensional analytical prediction of the dE/dx curve for different nuclear 
elements of the gas mixture is shown in colored lines. 
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Figure 12: Left: A fit to a recoil spectrum (markers) using a likelihood model based on 
helium and fluorine recoils. The helium contribution (dashed) and total recoil distribution 
(full line) are shown as histograms. Right: Unfolded neutron spectrum from a calibration run 
with 252 Cf source (markers), and a spectrum based on the 252 Cf table [29] (histogram). 



9. Future Prototype Underground Measurement 

We are proposing to use the prototype to make a measurement in the far-hall 
of Double Chooz at 300 m.w.e. We can use the 13 neutrons per day measured 
in the surface run to make an estimate for our event rate in the far-hall. For 
the surface run, we measure a muon rate of 100 muons per 2 s 1 . We assume 
the mean muon energy is ~3 GeV. The muon rate in the far-hall is 0.4 muons 
per m 2 s 1 with a mean muon energy of 60 GeV [4]. Using the power law scaling 
from Ref.[20 , we estimate ~0.5 events per day. 

A study of the prototype's efficiency as a function of the incoming neutron 



energy is shown in Fig. 13 In the study we assume that fluorine and helium 
recoils are efficiently reconstructed above 150 and 100 keV of recoil energy re- 
spectively. Since fluorine recoils are shorter, they have a better chance of fitting 
into CCD view-field. Geometric acceptance efficiencies are shown as full lines 



in Fig. 13 Helium recoils receive more energy from a collision with neutrons 
and leave longer tracks - a combination that is favored by reconstruction al- 
gorithms. Helium recoils are more efficient at energies < 20 MeV. Fluorine 
appears more efficient if the detector is being used for neutrons with the energy 
of tens of MeV's due to more compact tracks. This study indicates that we will 
be sensitive to neutrons from (a, n) and muon interactions described in Section 

o 



9.1. Conclusion 

We have modified the DMTPC dark matter detector design to produce 
a 2.8 L prototype directional neutron detector. The primary change to the 
DMTPC design was to add helium as a target for the neutrons. We have 
demonstrated that a gas mixture of 75 Torr of CF4 and 525 Torr of helium 
is capable of extracting the energy and angular distributions of the neutrons 
from a 252 Cf source. Since nuclear recoils are mostly produced at lower energies 
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Figure 13: The efficiency for reconstructing helium and fluorine recoils as a function of neutron 
energy. The geometric efficiency (full line) requires full track containment in the CCD view- 
field. The full efficiency (dashed line) assumes a step-function for the recoil reconstruction 
efficiency with thresholds given in the text. 

and have a higher probability of being contained in the CCD view-field at these 
lower energies, sensitivity to higher neutron energies is possible. 

The sensitivity to higher energy neutrons and an estimate of 0.5 events per 
day at the Double Chooz far site indicates that the prototype could make an 
interesting measurement at the far site on its own. The success of this prototype 
motivates the construction of a full-sized DCTPC system. This will consist of 
two 60 1 detectors, located at the near and far halls of Double Chooz, which 
can be interchanged for comparison of systematics. The goal is to measure the 
neutron flux, with information on both energy and direction, as a function of 
time, in order to tune the fast neutron simulation. 
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